EPJ manuscript No. 

(will be inserted by the editor) 



OO 
On 
On 



£s Low temperature thermodynamics of inverse square spin models 
< in one dimension 



On 



CO 

a 

C 

o 
o 



W. Hofstetter 1 and W. Zwerger 2 

1 Theoretische Physik III, Universitat Augsburg, 86135 Augsburg, Germany 

2 Sektion Physik, Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, 80333 Miinchen, Germany 

April 1998 

Abstract. We present a field-theoretic renormalization group calculation in two loop order for classical 
0(7V)-models with an inverse square interaction in the vicinity of their lower critical dimensionality one. The 
magnetic susceptibility at low temperatures is shown to diverge like T~ a exp(6/T) with a — (N — 2)/(N— 1) 
and b = 2n 2 /(N - 1). From a comparison with the exactly solvable Haldane-Shastry model we find that 
the same temperature dependence applies also to ferromagnetic quantum spin chains. 
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1 Introduction 

The investigation of long range forces in statistical me- 
chanics has a considerable history, including studies on 
the basic requirements for a proper thermodynamic limit 
jj| or on the existence of phase transitions even in one- 
dimensional systems. For classical spins, it has been known 
for some time that ferromagnetic models with a pair- 
interaction decaying like n,~( 1+CT ' (a > is necessary for 

| a proper thermodynamic limit) exhibit no transition if 
a > 1. Choosing a < 1 however, it turns out there is in- 
deed a phase transition even in one dimension both for 
Ising models || as well as for systems with a continuous 
symmetry [Q. The borderline case a = 1 of an inverse 
square interaction is thus of particular interest. For Ising 
spins an early conjecture by Thouless || that this model 
exhibits a peculiar continuous transition with a finite jump 
of the order parameter was finally proven rigorously by 
Frohlich and Spencer Q . In the case of a continuous sym- 
metry like the XY- or Heisenberg-models with inverse 

. square interaction, Simon has shown |jj that no symmetry 
breaking appears here, suggesting that the borderline in- 
teraction exhibiting a Kosterlitz-Thouless type transition 
is (hin)/n 2 . An alternative proof that the inverse square 
XK-model in one dimension exhibits no long range order 
was later given by Simanek || , using the classical version 
of the Mermin- Wagner theorem. 

A quantitative calculation of the low temperature sus- 
ceptibility of long range ferromagnetic spin models was 
given by Kosterlitz || within a momentum shell renor- 
malization group (RNG) up to one loop order. In partic- 
ular the susceptibility was found to diverge exponentially 
like exp(27r 2 / (N — 1)T) in the case a — 1 of an inverse 
square interaction for models with a continuous symme- 
try, where only the trivial fixed point at T = describing a 



fully ordered ground state exists. Independe ntly , in a brief 
note, Brezin, Zinn- Justin and Le Guillou ]10[ presented 
the results of a two loop field-theoretic RNG of the O(N)- 
model with interaction proportional to n~( 1+a \ They de- 
termined the critical exponents up to order (d — a) 2 in 
the vicinity of the lower critical dimension d c = a. Here 
we present a detailed version of the field-theoretic two 
loop calculation which essentially follows their ideas. How- 
ever, it uses a somewhat different method to determine 
the renormalization factors and, moreover, concentrates 
on the behaviour in the marginal case d = a = 1 . 

The study of one-dimensional classical spin models 
with long range interactions is interesting not only from a 
pure statistical mechanics point of view, but is also rel- 
evant for actual physical problems. In fact the inverse 
square Ising-model is closely related to the well-known 
Kondo-problem of a magnetic impurit y lo cally coupled to 
the conduction electron spin density pl| . A more recent 
application of the same model is the two state system with 
ohmic dissipation, arising e.g. in the problem of quantum 
coherence between macroscopically different states |l2],[l3| . 
For the inverse square XY-c&se it turns out that the prob- 
lem of strong tunneling in a metallic single-electron-box 
can be treated by calculating the free energy of the clas- 
sical spin model in the presence of twisted boundary con- 
ditions Q. It was the latter problem which motivated 
our present work, whose intention is in fact twofold: First 
of all, we give a detailed account of the two-loop RNG 
of the 0(N)-model with inverse square interaction which 
has obviously not been done in the literature so far. In 
addition to that we discuss the corresponding S — 1/2 
quantum spin chain. In particular, it is found that the low 
temperature behaviour of the susceptibility is essentially 
identical with that for classical spins. 



2 RNG in two-loop order 
a) Model and Definitions 

We start with a model of classical unit spins S(n) with 
N components which are placed on a ring with L sites, 
thus enforcing periodic boundary conditions. The inter- 
action energy in an external magnetic field H along the 
iV-direction is given by 

H cl = -J2[d(n-n')]- 2 S(n)S(n')-Hj2S N (n) (1) 



with 



din) 



L 7r|n| 
— sin 

7T L 



(2) 



the chord distance between spins which are n sites apart. 
The associated partition function can then be written as 

Z = IU n dS(n) 5{S 2 (n)-l)cxp{-S(T,H,S(n))}. (3) 

with an action 



S 



-E 



(S(n) - S(n')Y 



2T ^ {din - n')) 2 



H 



(4) 



Here J dS(n)S(S 2 (n) — l) denotes an integration over all 
orientations of a classical iV-component unit vector S(n). 
As is known from the work cited above, this model ex- 
hibits a continuous phase transition to an ordered phase 
below a critical temperature T c ~ t/{N — 1) for noninte- 
ger dimensionality d = 1 + e (e > 0). In order to study 
the behavior near the lower critical dimensionality one, we 
can therefore parametrize the spin S according to a low 
temperature expansion 



S(x) = n(x),v/i-n 2 ( a 



(5) 



The Hi (i — 1, . . . , N — 1) are Goldstone modes corre- 
sponding to the fluctuations around the ordered state with 
Sn = 1 where all spins point in the iV-direction. The gen- 
erating functional is now given by 



Z 



n 7= = ; exp{-5(T, H, H(n))}. (6) 

The term in the integration measure 

TT ; 1 : = exp|-lyin(l-n 2 (n))i (7) 

is due to the ((-function constraint and ensures the O(N)- 
symmetry of the model, which appears to be broken by the 
above parametrization. In a continuum notation Y] — * 
a -d j ( jd x w j^ n a lattice constant and 

TT . 1 ==expj-^a- d / d d x\n(l-U 2 (x))\ 
A „ A Jl-n 2 (n) l \ 2 J V y 



(8) 



In the following, we will work in the dimensional regular- 
ization scheme, where 



.-d 



d d q 
(2n) d 



(9) 



We can therefore neglect the measure term (see |iq| ). 

In order to generate correlation functions of the II 
fields, we introduce additional source fields according to 



S^S - dxJ(x) U(x). 



(10) 



The quantities which are most conveniently calculated di- 
agrammatically are the irreducible vertex functions I^™- 1 . 
Their generating functional r[TL, H] is formally obtained 
from lnZ[J, H] by a Legendre transformation with respect 
to J where II is the corresponding classical field. In our 
calculation, we will use the two point function 



r {2) {x iy ) 



s 2 r 



sn(x)sn( y ) 



(ii) 



As in the case of the non-linear cr-model with next- 
neighbour coupling jl5 16 it can be shown that the model 
is renormalizable by introducing two renormalization fac- 
tors Z t and rescaling the temperature and the field II. 
These factors must be adjusted in such a way that the 
renormalized iV-point functions 



(12) 



remain finite in terms of the renormalized parameters 

t = n e Z i T 1 T (13) 



h = Z^ZZH. 



(14) 

Below we will see that the choice Z t = Z v = Z with a 
single independent renormalization constant Z is possible. 



b) Details of the two-loop-calculation 

In order to calculate the Z-factors introduced above to the 
order i 2 , we first rewrite the action in momentum repre- 
sentation: 



S 



1 



d d q 



U 2 (q) (U 2 ) 2 (-q) 
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i%)n(- 



o((u 2 y 



U 2 (q)U 2 (-q) 



(15) 



where we have rescaled H i-> 2irH and T i-> {2n) d T. The 
free propagator corresponding to the quadratic term has 
the form 



(II a (p)II0(q)) =6 af3 5(p + q)G%(p) 
= (2ir) d 5 al3 6(p + q) 



(16) 
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Fig. 1. Diagrammatic representation of the four- and six-point 
interaction terms. 
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Fig. 2. Free propagator (a)) and one-loop contributions to 
T (2) (b) and c)). 



Fig. 3. Two- loop contributions to P™. 



with a, [3 = 1, . . . , N — 1. In addition, S contains two 
interaction terms which contribute up to two-loop order. 
These are shown in fig. ([!]) and correspond to a four- and a 
six-point-interaction. We will renormalize the model using 
the two-point function r^ 2 \ The relevant diagrams up to 
two-loop order are shown in fig. (||) and (||). In a short 
hand notation, the two point function in terms of these 
diagrams is given by 

r^(k) = a ) + b) + c)-d)-e)-f)- g) - h) - i) 

+j) + k) + I) + m) + n) + o) + p) + q). (17) 

Each one of these diagrams is proportional to T L ~ X where 
L is the number of loops. Several of them contain diver- 
gences for e — > 0. In the minimal subtraction scheme, we 
want to absorb these divergences in the two renormal- 

(2) 

ization constants Z t and Z v . After renormalization, 
should be an analytic function of the variables h, t and p. 

For the action (|l5|) we note that the p-dependent parts 
of the diagrams are finite. To one-loop order, for instance, 
only the second diagram (labeled by c)) depends on the 
external momentum 



d d q l + \p-q\ 



d d q \p - q\ 



{2n) d l + \q\ J {2n) d 1 + \q\ 
Evidently this integral is finite for d < 2. Because of 



(18) 



4 a) (p) = i L « 
At 



e \p\ + Z t Zn 



+ higher order terms, 



(19) 

the coefficient of \p\ in the inverse propagator is renormal- 
ized by Z^/Zt- As is evident from the one- loop term (|l8|), 
this prefactor is already finite. Therefore the choice 



Z JY — Zi 



z 



(20) 



is possible, leaving only one independent renormalization 
constant. In fact this relation is valid to arbitrary loop or- 
der, as was shown by Brezin, Zinn- Justin and Le Guillou 
]l0| ]. As a consequence, we can take the external momen- 
tum p to be zero in the following. One further simplifica- 
tion is due to the fact that in dimensional regularization 



"empty" integrals of the type J d d q are equal to zero. This 
implies that the contributions of the diagrams in c), h), n) 
and p) vanish. Moreover, several contributions turn out to 
be identical and therefore cancel: this applies to the dia- 
grams d) and o) which are given by 



TH 



2d-l 



d d qd d r l + \q 



(21) 



(2ny d (i + MXi + M)- 

Similarly the diagrams i) and m) which give a contribution 



AT- 1, 



-TH 



2d-l 



d d qd d r 
(2ir) 2d 



1 



(l + M)(l + |r|) 



(22) 



are equal and therefore cancel in r^ 2 \ The remaining ex- 
pression to be evaluated is 

(o, H) = a) + b) - e) - f) - g) +k)+j) + l)+q). (23) 

In order to extract the singular parts of the integrals, we 
use the asymptotic expansions 



1 



and 



d d q 
{2ir) d l + \q 

d d q 



K d 



1 



+ o(e) 



1 



(2tt)<* (1 + M) 2 
where the constant 



= KA\ 



2 )) 



(4n) d / 2 r(d/2) 



(24) 



(25) 



(26) 



results from the d-dimensional angular integration. The 
identities (|24|) and 
gularities in e for al 



allow to extract the leading sin- 
contributions to d23| ) except e). By 
a detailed asymptotic analysis of the corresponding inte- 
gral, the following singular contributions of this diagram 
can be obtained: 



e) 



d d qd d r 



1 



<1\ 



(27T) 



2d 



N - 1 



TH 1+2t K 2 , 



?l) 2 (l 



M) 
i 

27 



(27) 



Neglecting nonsingular parts, the resulting two-point func- 
tion is therefore finally given by 



T 2e 



TH 



l + 2c 



(N — l)(JV-2) 



4c 



3(JV- l) 2 



(28) 



where we have absorbed the angular integration term Kd 
in the temperature according to 



T^TK 



r 



(2) 



(29) 



At this stage we can now introduce renormalized parame- 
ters h and t as defined in (JL3|) and (|l4|). Taking k = 1 for 
convenience, we obtain 



4 2) (M) 



zn N-i 



+th 



— h 1+ 

t 2e 

l+2e (N-l)(N-2) , 



th 



(30) 

I 1 2e 3(iV-l) 2 



4e 

where terms of order t 2 have been neglected. Using the 
ansatz 

Z= I + ai (e)t + a 2 (e)t 2 + o(t 3 ) (31) 
this can be written as 

'4a 2 - a\ 



— \- n \ — ai 
t V2 



TV — 1 



2e 7 V 8 
(TV - 1)((3 + e)ai - (N — 2)) 3(7V - l) 2 



+i In ft 



4e 8e 2 
3(7V-1) 2 (iV-l)(3 + e) 
4e 2 



-ai 



(32) 



where we have neglected terms which are not singular in 
e. The coefficients a\ and a 2 are now chosen such as to 
eliminate the poles of order e -1 and e~ 2 . This is achieved 

by 

oi = ^— ^ (33) 



and 



a 2 



(N — l) 2 (JV-1) 
e 2 2e : 



(34) 



thus rendering the renormalized two point function finite. 
The renormalization constant Z is therefore given by 



Z = 1- 



iV-1 /(iV-1) 2 jV-1 



2c 



t 2 + 0(i 3 ). (35) 



From this result we can now determine the beta function 

(36) 



measuring the variation of the renormalized coupling un- 
der variation of the momentum scale, keeping the bare pa- 
rameters T and H fixed. According to (|l3|), we can write 



/?(<) = -et ( 1 + t 



.dhxZ 
dt 



(37) 



Using the result ([35]) this immediately leads to 

(3{t) = -et + (N — l)t 2 + (N - l)t 3 + 0(t 4 ). (38) 

For e > 0, the renormalization flow has a non trivial fixed 
point corresponding to a phase transition at finite tem- 
perature. The associated critical behaviour has been dis- 
cussed by Brezin et al. |l(J . Here we are interested in the 
one-dimensional case, where only the trivial fixed point 
at t = describing the fully ordered ground state ex- 
ists. During the calculation we have rescaled the coupling 
T twice. Inserting the original temperature according to 
T — > T/2ir 2 we get the following equation for the flow of 
T under the reduction of the cutoff A — ► Ae~ l : 



dT 



rp2 rp3 



0(T 4 ). (39) 



Integrating this equation between an initial temperature 
T very close to zero and a final value T = O(l), the cor- 
responding length rescaling factor behaves like 



£(T) - T~ exp 



2tt 2 



(N - 1)T 



(40) 



Obviously £(T) plays the role of a "correlation length", 
however it is important to point out that there can be 
no exponential decay in models with long range interac- 
tions. In fact, using Lieb-Simon-type inequalities, it can 
be shown that the spin-spin correlation function 



g(n) = (S(n) S(0)) 



(41) 



asymp totically decays like n~ 2 for arbitrary temperatures 
|L7jM - In spite of the absence of an exponential decay, 
within a standard scaling hypothesis the behavior of g(n) 
near the critical point at T — is described by a universal 
function f(n) with /(0) = 1 and f(\n\ — > oo) = f 00 /n 2 1 
foe — 0(1), such that 



g(n,T^0)=f 



( n 



(42) 



with £(T) given by (|40|). As a result the magnetic suscep- 
tibility 



X 



= /3 l + 2]Tff(7 



(43) 



is finite for any T and scales like 



X(T) ~ T~— exp 



2tt 2 



(N — 1)T 



(44) 



at low temperatures. The two loop calculation thus 
uniquely determines both the exponent and the In- 
dependence of the prefactor of the magnetic susceptibility 
as T -> 0. 
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3 Discussion 

In this work we have calculated the low temperature sus- 
ceptibility of classical O(N) spin models with an inverse 
square interaction. Now it is interesting to compare our re- 
sults with the corresponding ferromagnetic quantum spin 
models. In the case of a Heisenberg model 



H = - J (n ~ n') S„S„, - H £ S z n 



(45) 



with interaction J(n) = d~ 2 (n), the exact low tempera- 
ture susceptibility in the case S — 1/2 was calculated by 
Haldanc p9|. It is given by 



X(T) 



4^ 



exp 



7T" 

4T 



(46) 



and thus exhibits precisely the same temperature depen- 
dence as our result ([id]) for the classical case with N = 3 
(note that we have chosen units in which J — I and our 
temperature differs from the one in ref. jl9| by a trivial 
factor two arising from ^2 n ^ n , — 2^ n <n')- ^ n ^ ac ^ the 
only difference is in the numerical factor which appears in 
the exponent. On a qualitative level the exponential diver- 
gence of x{T) in the case of long range interactions can be 
understood from a simple spin wave calculation around 
the perfectly ordered ferromagnetic ground state. For a 
general spin S, the Hamiltonian ([l5|) leads to a dispersion 
relation of magnons which has the form 



e(k) = 4S ^ J(n) (1 - cos(fcn)) + H. 



(47) 



n>0 



An equilibrium distribution of noninteracting magnons 
leads to the standard reduction in the dimensionless mag- 
netization 



m(T,h) = 1 



elk 



1 



2tt e^W - 1 ' 



(48) 



For the Haldanc-Shastry model with J(n) = d 2 (n), the 
magnon dispersion (|47|) is given by 



e(k) = S(2n\k\ - k 2 ) + H 



(49) 



(see also 19 ). The associated reduction in magnetization 
in the limit H, T — » then behaves like 



Sm(T, H) 



T 



2(ttS) 2 



hi 



2tt5 
H ' 



(50) 



The condition 8m — b with a constant b of order one then 
defines a field H(T) which qualitatively determines the 
zero field susceptibility via 



X(T) 



1 



H(T) 



exp • 



2b(irS) 2 



T 



(51) 



Choosing b = 1/2, this result agrees precisely with the 
exact solution (|46| ) for S = 1/2, while the exponent in the 



susceptibility of the classical spin model with unit spin is 
obtained by formally setting S = 1 . A modified spin wave 
theory and Schwinger boson techniques have been applied 
to the Haldane-Shastry model by Nakano and Takahashi 

ill- 

For short range interactions where J(n) decays faster 
than n -3 , the magnon dispersion has the usual form 



e(k ->0) = Dk 2 



H. 



(52) 



It is then straightforward to see that the above spin 
wave calculation leads to a susceptibility diverging like 
%(T) ~ S 2 /T 2 . This behavior is indeed found from the 
exact solution of both the nearest neighbor S — 1/2 case 
[p2| as well as for the corresponding classical Heisenberg 
model Q . We thus expect that for any ferromagnetic spin 
model with a fully ordered ground state the temperature 
dependence of the susceptibility at low T is independent 
of whether the spins are classical or quantum mechanical. 

In the case of an AF-modcl (N = 2), the ground state 
for the S = 1/2 chain with inverse square interaction has 
again a very simple structure, as shown by Haldane p4[ . It 
would thus be interesting to check whether in this case the 
low temperature susceptibility again diverges like that of 
the classical model with N = 2, i.e. x{ T ) ~ exp(27r 2 /T). 
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